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ABSTRACT 

A general purpose coiaputer program for the 
calculation of a matrix of tetracho::ic correlations is described. 
This program was developed for use in adaptive (ana other) testing 
research for examining the unidimensicnality assumption in latent 
trait theory, in conjunction with available tactor analysis programs. 
Several other potential applications, as well as details for its use, 
are described. The program accepts as input raw dichotomous data, 
reduced joint frequency data, or joint and marginal proportions, for 
up to 75 items. Output options include the tetrachoric correlation 
matrix^ the matrix of phi coefficients, fourfold frequency tables for 
every item pair? a joint frequency matrix (which reduces all the 
information in the fourfold tables to a square matrix with order 
•equal to the number of items) ^ and a pair-by-pair listing "^f input 
proportions and output correlations which permits testing the program 
against published tables of the tetrachoric correlation. Variable 
input anc** output formatting makes the program convenient to use in 
conjunctions with other analyses by packaged statistical programs. 
Examples of input and output are presented. A complete FORTRAN IV 
listing is included. (Author) 
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Disclaimer 



The computer program presented In this report has been 
tested carefully In the analysis of a variety of data. Although 
It has consistently yielded accurate results, the authors make 
no warranty as to Its accuracy and functioning, nor shall the 
fact of Its distribution Imply such warranty. 



TETREST: A FORTRAN IV PROGRAM FOR CALCULATING 
TETRACHURIC CORRELATIONS 



The Tetrachorlc Correlation Coefficient 

The tetrachorlc correlation coefficient Is an estimate of the product- 
moment correlation between two continuous normally distributed variables which 
have been dichotomized so that each variable takes only two discrete values. 
Stated another way, the ^.etrachorlc r Is an estimate from dlchotomous data of 
the blvarlate normal correlation between two continuous variables which have 
been dichotomized. Use of the tetrachorlc correlation Involves assumptions 
that 1) a continuous latent variable underlies the dlchotomous data on each 
variable; 2) the latent variables are normally distributed; and 3) the 
regression of one on the other Is linear. 

These assumptions are Illustrated In Figure 1, whl^:h deplrts normally 
distributed variables and which have been dichotomized at points g^^ 

and g^t respectively. In such a way that If the continuous value of X^ equals 

or exceeds Q • a binary value of 1 Is assigned; if X Is less than g a value 

n n 

of 0 is assigned. Lines drawn perpendicular to the axes of X^ and X ^ through 



points g^ and g^p, Intersect to divide the blvarlate normal scatterplot of 



Figure 1 



Partloning of a Blvarlate Normal Distribution Scatterplot 
into Four Quadrants by Dichotomlzation 
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and X2 into four quacirants labeled Cy d. 

The tetrachoric correlation is used to infer the bivariate normal 
correlation from the frequencies observed in each of the four quadrants. 
These frequencies are the joint 'frequencies wliich constitute a fourfold, or 
2x2 contingency table of the two variables. Frequency a is the number of 
times that a value of 1 was observed on both variables jointly; b is the 
number of occurrences of a 1 on jointly with a 0 on X^; c is the frequency 

of a (0,1) pattern on X^^ and X2, respectively; d is the frequency of occurrence 

of the (0,0) pattern. Figure 2 illustrates the correspondence between the 
partitioned bivariate scatterplot and the fourfold table; each of the frequencies 
Uy by Cy d corresponds to the observed frequency in the same-lettered quadrant 
of Figure 1. 

Figure 2 

The Fourfold (2x2) Table of Frequencies 
Corresponding to the Quadrants in Figure 1 




a+b 



c+d 



b+d a+c 

Figure 2 also depicts the marginal, or unconditional, frequencies of the 
two variables. For example, a-fb is the unconditional frequency with which a 
value of 1 was observed on variable X^^; ahe is the unconditional frequency of 

a 1 on X^. The marginal frequencies of a variable n are a function of the 

point of dichotomy, g . If the assumptions underlying it are met, the 

n 

population tetrachoric correlation coefficient is invariant under changes in 
the point of dichotomy. This is not true of the phi coefficient, the only 
other correlation index used widely with dichotomous test items. The magnitude 
of phi is heavily depetident on the relative magnitudes of the marginal frequencies 
of the two variables. 

The known lack of invariance of phi when there is diversity in the marginal 
frequencies from one item to another has led to increased use of the tetrachoric 
correlation. Whether tetrachoric correlation coefficients are Invariant 
regardless of the points of dichotomy, when the assumption of bivariate 
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normality is not met> is an empirical question whose answer is specific to the 
problem at hand. For bivariate normal data, the tetrachoric correlation 
estimates the population product-moment correlation and is invariant. For 
data which do not meet the assumption of bivariate normality the tetrachoric 
correlation can still be computed, but should be interpreted only as a measure 
of association whose invariance under disparate marginal frequencies is 
questionable* 

Phi coefficients and tetrachoric correlations also differ in another 
important respect. A matrix of sample phi coefficients is always non-negative 
definite, or Gramian. This means that such a matrix is appropriately structured 
for the application of certain statistical analysis techniques, including 
factor analysis. Matrices of sample tetrachorics are often non-Gramian 
(Lord and Novick, 1968, p. 349). Factor analysis of a non-Gramian correlation 
matrix of full rank will result in one or more factors having negative eigenvalues; 
that is, "factors" which account for negative variance. Related problems may 
include communality estimates exceeding unity, factor loadings exceeding unity, 
and inability to estimate communalities iteratively because the inverse of the 
correlation matrix is indeterminant . 

The non-Gramian property of some sample tetrachoric correlation matrices 
may be due to violating the normality assumptions for one or more variables, 
to sampling error where the population interitem distributions are all bivariate 
normal, or to numerical errors in estimating the correlations. Any or all of 
these sources of error may render hazardous the interpretation of factor 
analyses based on tetrachoric correlation matrices. Users of the tetrachoric 
correlation should be aware of these potential problems, and should exercise 
caution in statistical analyses of tetrachoric correlation matrices. 

Applications of the Tetrachoric Correlation Coefficient 

Most recent research in adaptive testing has employed a latent trait 
theoretical model to account for testees* responses to dichotomous test items. 
The most prevalent latent trait models assume that the latent variable under- 
lying test item responses is one-dimensional. For constructing test item 
pools, this assumption requires that the latent space spanned by the items be 
unxdimensional. A sufficient condition for the unidimensionality of the latent 
space is that the matrix of tetrachoric item intercorrelations has just one 
common factor (Lord and Novick, 1968, pp. 381-382). For practical purposes, 
the items are considered to represent a unidimensional variable if the largest 
latent root of the sample matrix of their tetrachoric intercorrelations 
accounts for a large proportion of the common variance, and the second and 
smaller latent roots are of about the magnitude that might be expected from 
sampling errors (Indow & Samejima, 1966; McBride & Weiss, 1974). 

The tetrachoric correlation, then, is a useful statistic in the develop- 
ment and analysis of item pools for adaptive testing. But it is useful for 
other applications where dichotomously scored items are used, by virtue of 
its invariance properties. If the assumptions underlying it are met, the 
tetrachoric correlation coefficient is invariant under differences in the 
marginal proportions of the two variables involved. As discussed above, this 
invariance property is not the case with the fourfold point (phi) correlation 
coefficient, which reaches its maximum value only if both variables of a pair 
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have the same marginal proportions (Carroll, 1961). For this reason, the 
tetrachoric correlation is useful in factor analyses of dichotomous variables 
whose marginal proportions vary widely. 

Another application of the tetrachoric v is in the estimation of the 
parameters of the item characteristic curve in latent trait t<»st theory. 
If all the items in a test are considered to have item characteristic curves 
of normal ogive form (Lord & Novick, 1968, Chapter 16), Bock and Wood (1971) 
have pointed out that under certain conditions each item's discrimination 
parameter may be estimated from its loading on the first factor extracted from 
the matrix of tetrachoric interitem correlations. 

Despite the increase in applications of the tetrachoric correlation, 
general purpose computer programs for calculating it are not widely available. 
The computer program described in this report was written specifically to fill 
this void. Its original application permitted a unidimensionality analysis 
of a set of test items being assembled in order to implement an adaptive 
ability testing program (McBride & Weiss, 1974). It is presented here in 
order to make a tetrachoric correlation program available to users who would 
otheirwise not have access to one. 

General Description of TETREST 

This program (TETRachoric ESTlmation) uses an approximation procedure 
given by Kirk (1973) to estimate the tetrachoric correlation coefficient 
between two dichotomized or dichotomous variables. It also calculates the 
value of the fourfold point (phi) correlation between two dichotomous variables. 

TETREST is designed primarily to: 

1. Construct a set of four fourfold (2 x 2) contingency tables, for 
every pair of variables, from dichotomous data on a set of K variables. 
For K variables there are [(K)(K^l)]/2 such variable pairs. 

2. Estimate the degree of relationship (tetrachoric and/or phi correlation) 
in each fourfold table. 

3. Construct the KajK matrix of tetrachoric and/or phi correlation 
coefficients in a form amenable to factor analysis using packaged 
computer statistical routines such as SPSS: Statistical Package 

for the Social Sciences (Nie, Bent and Hull, 1970) or HMD: Biomedical 
Computer Programs (Dixon, 1973). 

The program is general, emphasizing flexibility and ease of input and 
output for the user. It was written specifically for the CDC CYBER 74 computer 
in Control Data FORTRAN IV. TETREST will accept input data from punched 
cards, magnetic tape, or disc storage. It will print, punch, or write its 
output on logical units in several different forms. The specific input and 
output options are described below. 

Data Input 

TETREST will accept input data of three different types, all arising from 
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dichotomous data on K variables observed over N individuals. The three data 
types are illustrated in Figures 3, 4, and 5 respectively. 

Type I is the basic form of data: the dichotomous scores of N individuals 
on K items. Figure 3a illustrates the Type I data of 10 individuals on three 

Figure 3 

(a) . Type I data from a 3-item test for 10 individuals 



Individual 


1 


Item 
2 3 


1 


1 


0 


0 


2 


1 


1 


0 


3 


1 


1 


0 


4 


0 


0 


1 


5 


1 


i 


0 


6 


0 


0 


1 


7 


1 


1 


0 


8 


1 


0 


1 


9 


1 


0 


1 


10 


0 


0 


1 


Totals 


7 


4 


5 



(b). 



The fo urfold (2x2) contingency tables for each 
pair of variables, derived from Type I data 



Item 1 



1 
0 



Total 



Item 2 
0 1 



Total 



10 



Item 1 



1 
0 



Total 



Item 3 
0 1 



5 
0 



Total 



10 



Item 2 



1 
0_ 
Total 



Item 3 
0 1 



Total 



4 
6 



10 



dichotomous items (variables). Each individual's item scores (1 or 0) are on 
a separate data record, such as a punched data card. Counting the frequency 
of I's and O's permits construction of the fourfold tables, shown in Figure 3b. 

Type II data is the reduction of Type I data to a KxK matrix of frequency 
counts. Figure 4 shows the Type I data of Figure 3 reduced to a 3 x 3 matriA, 
called the joint frequency matrix, which contains the minimum frequency 
information needed to construct the fourfold tables for each of the [(K)(K-l)]/2 
possible pairs of variables. The three diagonal entries (7,4,5) are simply 
the three marginal "positive" frequencies (e.g., the number of correct or 
endorsed items). The supra-diagonal entries (4,2,0) are the joint positive 
(i.e., 1,1) frequencies for variable pairs (1,2) (1,3) and (2,3) respectively; 
that is, the frequency in cell a of the fourfold table of variables i and J. 
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The infra-diagonal entries (3,0,1) are the joint negative frequencies 
(i.e., 0,0) for variable pairs (2,1), (3,1) and (3,2) respectively; that is, 
the frequency in cell d of the fourfold tables (see Figure 2). 



Figure 4 



The Joint Frequency Matrix of Type II Data 



1 

Item 2 
3 



7 4 
3 4 
0 1 



Figure 5 shows Type III data redu'-ed from the fourfold tables in 
Figure 3b. Type III data consists of three proportions (P... a.* a J and 

and optional sign-change character, SIGN. Each set of three proportions is 
derived from the fourfold table of one variable pair. 

is the joint proportion obtained by dividing one of the cell frequencies 



(a^hyC ov d) by Ny the number of individuals observed (N=a'fb'fa'i'd) . The 
appropriate cell is determined by the marginal proportions. P^^ is the joint 

proportion occurring in the cell corresponding to the smaller marginals. 

q. is the smaller marginal proportion for variable 

q. is the smaller marginal proportion for variable j. 

J 

SIGN is a character which changes the sign of the obtained correlation. 
This is necessary whenever P'j comes from cell b ox c of the fourfold table. 

If SIGN is left blank, no sign change will occur. Any other character, such 
as a minus sign (-) , will effect the sign correction. 

Figure 5 

The Joint Proportions Matrix of Type III Data 



Variable 
Pair 
(i.j) 




^i 




SIGN 


(2.1) 


0 


.40 


.30 




(3,1) 


.30 


.50 


.30 




(3.2) 


0 


.50 


.40 





Data of Type I will likely be used most often in calculating tetrachoric 
correlation coefficients. TETBEST* is capable of yielding data of Tyre II as 
output since some program users might desire to use such output as input at a 
later time. ^?here data of Type II are available, their use as input can 
represent a considerable saving of computer processing time, especially when 
the number o'c individuals in the original data is very large. The provision 
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for data of Type III permits rapid calculation of tetrachoric or phi correlations 
for one or more variable pairs for which raw data or frequency matrix data 
are not available. One application of Type III data might be to check TETREST 
results against tables of the tetrachoric correlation coefficient, such as 
those published by the National Bureau of Standards (1959). 



Data Output 

TETREST is capable of several output options; the output ? \i depends 
in part on the type of input data used. 



For Type I and Type II data, the following output options are available: 

1. Printing or writing on tape or disc of the KxK joint frequency matrix 
(i.e.. Type II input). 

2. Printing of the [(K) (K''l)]/2 2x2 frequency tables, including values 
of phi and tetrachoric v. 



3. Printing, punching, or writing on tape or disc or the KxK matrix of 
tetrachoric and/or phi intercorrelations in either lower diagonal 
or square format. 



For Type III data, option 3 is available. Also available optionally for 
Type III data is a listing by variable pair of the input proportions, tetra- 
choric and phi correlation values, and the degree of convergence . occurring 
in the approximation of the tetrachoric r« 



Program Use 



TETREST incorporates a number of input and output options which may be 
specified by the user. Implementing these options involves the use of from 
five to ten program control cards. Each card has a task definition field 
consisting of the first 15 columns of the card, and one or more specification 
fields spanning the remaining 65 columns (i&"80) . Program users familiar 
with the SPSS system (Nie, Bent, and Hull, 1970) will recognize that the 
present system of control cards is very similar to the one used within SPSS- 
The task definition field of each control card defines the nature of the 
data appearing in the specification fields. The ten program control cards* 
in the order in which they should appear, are as follows: 



PROBLEM 
RUN NAME 
INPUT MEDIUM 
INPUT FORMAT 
N0# OF CASES 
OUTPUT FORMAT 
OPTIONS 

READ INPUT DATA 

MATRIX 

FINISH 



(mandatory) 

(optional) 

(optional) 

(mandatory) 

(optional) 

(optional) 

(optional) 

(optional) 

(mandatory) 

(optional) 



Descriptions of the control functions of each card appear in sections A 
through J below. 
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A, The PROBLEM Card Cmandatory) 

This card may contain from 1 to 6 specification fields, as shown below: 
PROBLEM HVAR^-- DftTA«>> 2X2^- NrHT=>-- KEY«=- R> 

The data in each field of the PROBLEM card must appear in exactly the columns 
specified. 

NVAR=»_ jspecif ies the number of variables in the analysis. NVAR must be a 
right-justified integer number greater than 1 and less than 76. 
NVAR is the only mandatory specification on the PROBLEM card. 

DATA* specifies the type of input data and some output information. 

If DATA is not specified, Type I data is assumed. DATA may take any 
integer value from 1 to 5, or may be left blank. These specification 
values result in the following input /output options: 

DATA"«1 or blank: the program will read Type I dai.a. 

2 : the program will read Type I data and output the joint 

frequency matrix calculated. 

3 : the program will read Type II data. 

4 : the program will read Type II data and output the joint 

frequency matrix. 

5 : the program will read Type III data. 

2x2»9L causes the program to print the [(K) (K''l)]/2 2x2 frequency tables, 
one for every variable pair, only if Type I or Type II data was 
specified. Any other numerical value, or a blank, suppresses printing 
of the 2x2 tables. Included in each 2x2 table are the intervariable 
phi and tetrachoric correlation coefficients. (CAUTION: each 2x2 
table requires 6 lines of print. A 75-variable problem will print 
[(75) (74)]/2»2775 tables, which requires over 200 pages of print for 
the 2x2 tables alone.) 

NFMT« ^specifies the number of cards used to specify the input format. 
NFMT may be a right-justified integer number from 1 to 10. If NMT 
is not specified; it is set equal to 01, and the variable input format 
must be no longer than 65 characters. (See INPUT FORMAT below.) 

KEY=» ^specifies the data character which vill be considered the "correct" 

or "positive" value in the analysis of Type I data. For example, in 
scored test data the character 1 might be assigned a correct response; 
in questionnaire data a I might represent a "yes" response. KEY may 
be any alphanumeric character. Once KEY is specified, any input data 
character other than KEY will be considered as "incorrect" or "negative." 
For example, in (binary) data a 1 would be considered correct if 
so specified; and a 0^ blank, or any other character would be considered 
incorrect . 

If KEY is not specified, it is automatically set to 1 (a Hollerith 
constant). If data other than Type I data are used, KEY is ignored. 
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R* ^specifies the correlation coef ficient (s) to be used in constructing 

the output correlation matrix or matrices. Either tetrachoric or phi 
correlation matrices, or both, may be specified by a left-justified 
alphabetical entry to the right of R=», 

R*PHI specifies that only a matrix of phi coefficients will be calculated, 

R»TET or a blank field specifies that only a matrix of tetrachoric 
correlations will be calculated « 

R"«BOTH specifies that both matrices will be calculated. 

B. The RUH NAME Card (optional) 

This card allows the program user to print any arbitrary 65-character 
label at the top of the first printed page of the data analysis • The label 
must be punched in columns 15-80 of the RUN NAME card. The RUN NAME card, 
which takes the form shown below, is completely optional. 

RUN HftME ftHY ftRBlTRflRY SET DF CHftRftCTEf S 

h ? 3 4l5 • 7 tit toil tt|t3M t5»tn »1i W»nMla»27?t0tg3T3;U3 34 353tin3tM4<MU?43 44k5W47 4lM»5^ 

C. The INPUT MEDIUM Card (optional) 

This card allows the program user to input data from cards, magnetic 
tape, or disc storage file. The specification must begin in column 16 of the 
card, and must be one .of the three values shown below in the example. If an 
INPUT MEDIUM card is not used, the program will read the input data from the 
INPUT file, which ordinarily implies punched cards. 
INPUT MEBIUH CARD i-ead data from file INPUT 

DISC ) read data from file DATA 
TAPE / 

If INPUT MEDIUM specifies DISC or TAPE, the user must supply the input data on 
a file named DATA, and must position that file at the start of the first 
record of the data to be read. (The program does not rewind the file DATA. 
This is to permit separate analyses of sequential blocks of data on the same 
file if desired by user.) 

D. The INPUT FORMAT Card (mandatory) 

The program reads data in an Input format specified on this card by the 
user. The program uses a variable format, the form •'^f which varies with the 
type of data to be analyzed. The INPUT FORMAT card is mandatory for the first 
problem to be analyzed, but is optional thereafter unless the format changes. 

Type I data must be read in Al format fields, since each variable score 
will be compared with the one-character alphanumeric value of the KEY 
specification. 

Type II and III data (except SIGN) must be read in floating-point numeric 
fields (F-fields), since the data are numeric and will be read into floating- 
point arrays. SIGN must be read in an Al format. 
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The INPUT FORMAT specification begins in column 16 with a left parenthesis, 
(, and the input format, which begins on the INPUT FORMAT card and continues 
tnrough NFMT cards, as specified on the PI.IOBLEM card. The input format may 
extend to as many as 10 cards in length, inclusive of the INPUT FORMAT card. 
Continuation cards of the INPUT FORMAT may use the entire 80 card columns. 
The format specification must end with a right parenthesis, ). 

INPUT FORMftT <...,,.,..,) 

h ? 3 4h I 7 th 10 M M IS <I|17 It tt;tt7171?3?4|2^?s;7?t<?il0 3'PiM 34 3SMt3/3l3l40M^ 

E. The NO. OF CASES Card (optional) 

This card specifies the number of individual records of NVAR variable 
scores for the program to read in using, the INPUT FORMAT specified by the user, 
in the case of Type I data only. The specification field begins in column 16 
of the control card, and must be a right-justified integer number. E.g., 

m. or CASES 0192 

1 1 > 3 4 15 A-7_ilt 10 II ttlq M H Wtn H 1i 77 ?3 ?40 » 27 30 31 32^ 14 35 Jt^? Uji 40j4l 43 <3 *4i4S 41 47 4#<l» SlM^ M M Vl^J « Miojn t2 t3 Mjtt M 17 l^t TO n T3|n ^4 n Tt\n 71 H M| 

The program interprets this as 192 .subjects • The program will read exactly 
the number of cases ol Type I data specified, unless it first encounters an 
end-of-file on the input data file. In that case it resets the number of 
cases to the number of individual data sets read before the end-of-fiie was 
encountered. A message to this effect appears on the control-card page if an 
end-of-file is encountered. Program users may employ this feature to their 
advantage when estimating the number of cases if the exact number is not 
known, simply by using any reasonable over- estimate of the number of cases. 

If a NO. OF CASES card is not read, the number is automatically set to 
9999, and the program will read 9999 individual data records, or until it 
finds an end-of-file mark as described above. CAUTION: In Type I data, the 
number of cases is reset to 9999 at the end of each problem, and must be 
respecified using the control card in each succeeding problem when appropriate. 
The number of cases has no bearing on the solution of Type III data problems, 
since data are input as proportions. 

F> Ttie OUTPUT FORMAT Card (optional) 

The format in which the matrices of tetrachoric and/or phi correlation 
coefficients will be printed and/or written on the logical output unit may 
be specified using this card. The card has the form: 

OU TPUT FDRHAT < > 

f r? 3 4 |5 I 7 t h to n tth3 14 tS|l? \t\%7m21 23 " ?tUl 30 31 3?t33 M 35 31^7 31 39 40M> 4? 43 44i4S 46 47 4»l 50 V 52153 S4 SS ^hl H M MKl » U Mi« H ^0 71 I2\n M Ij ttVl H 

The specified format must be enclosed in parentheses, and must be contained 
entirely in columns 16 through 80 on the OUTPUT FORMAT card. The format 
should consist of F-fields, in order to provide for the floating-point numerical 
values of the correlation coefficients . 

If no OUTPUT FORMAT card is read, the program will write and/or print 
matrices in 8F10.7 format. 
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G. The OPTIONS Card (optional) 

The specification on the OPTIONS card determines the disposition of the 
correlation matrices resulting from the program* s calculations. The value 
specified must be an integer number between 1 and 6, inclusive. Theise values 
cause the following dispositions of the r-matrices: 

DPT I DNS 1 PUNCH the r-matrix on data cards 

I I ? 3 4 IS t 7 > 19 10 II 12in 14 >S I8|l7 H H 20 

2 PUNCH and PRINT the r-matrix 

3 PRINT the r-matrix on the output file 

4 PRINT the r-matrix and WRITE it on TAPES 

5 WRITE the r-matrix on file TAPB8 

6 suppress any output of the r-maOrix 



Defult value of the OPTIONS card is 3; i.e., pmnted output only. 

NOTE: 1. The OUTPUT FORMAT, whether set by the program or specified, 
applies only to printing or writing the ^r'-matrix. If PUNCH is specified 
(option 1 or 2), the PUNCH format is 8F10.7 whenever square matrix output is 
chosen (i^- iceeping with the SPSS factor analysis program's requirements for 
a square correlation matrix in that format) . When the lower diagonal matrix is 
selected (via the STATISTICS card), punched output is in format (10F8.5), 
which is somewhat more economical with little loss of accuracy. See section I, 
below, for infomiation on the square and lower diagonal matrix output options, 

2. For options 4 and 5, each correlation matrix is written as a file on 
TAPES, separated from other files by an end-of-file mark. When both tetrachoric 
and phi matrices are specified (via R«BOTH on the PROBLEM card) the tetrachoric 
matrix is written first, the phi matrix second. 

H. Tlie READ INPUT DATA Card (optional) 

This card has no specification field, and has the following form: 

READ INPUT DAT ft 

I I 2 3 4 Is < I h 13 n 14 tS]l7 H W ?o|21?? n ?4t?5 H V ?ii29 30 t 33|j3 M K Xjj? 31 38 40j4L4l 43 44l<S *i 47 4||4»05t 52|S3 S4 1,5 itj^l SI SI tO[%\ H U Ui<S if ^7 tifeljO n '2\n 74 n\l1 f> 19 toj 



Its function is twofold: 

1. It causes the program's internal data file (TAPEl) to be rewound 
(i.e., repositioned at the start of the first logical record), 
effectively erasing the input data from the previous problem, if any. 

2. It causes the program to read from the user's data file (INPUT or 
DATA) the input data of the uype specified on the PROBLEM card. 

For Type I data, the READ INPUT DATA card causes the program to read 
N (=N0. OF CASES) individual records of K (=NVAR) variable scores from the 
input data source specified, according to the INPUT FORMAT specified. If an 
end-of-f ile mark is read before the Nth data record , N is reset as described 
in section E, above. 

For Type II data, the READ INPUT DATA card causes the program to read 
O K («NVAR) records of K frequencies each (i.e., a KxK joint frequency matrix) 
according to NVAR as specified on the PROBLEM card, and in the INPUT FORMAT 

^ 17. 
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specified. If an end-of-file mark is encountered before the Kth full record 
is read, processing is terminated, and an error diagnostic is printed. 

For Type III data, the READ INPUT DATA card causes the program to read 
[(K) (K-DY/S separate sets of proportions, one set for each pair of variables, 
and later to process them as though they were infra-diagonal entries in a 
KxK matrix of proportion sets. Thus the first set results in the obtained 
correlation being treated as though it were between variables 2 and 1; the 
second between 3 and 1; the third between 3 and 2; the fourth between 4 and 1, 
and so on. The example below is for a 4-variable problem with [(4)(3)]/2=6 
variable pairs. 





Data 


Variable Pair 


1 


P q, q, sign 


(2,1) 


2 


P q, q, sign 


(3.1) 


3 


P qi q, sign 


(3,2) 


4 


P q, q, sign 


(4,1) 


5 


P q, q, sign 

t ^9 


(4,2) 


6 


p q. q sign 


(4,3) 



r-matrix: 




The numbers within the r-matrix schematic above correspond to the input data 
sets of proportions (see Figure 5). 

If no r-matrix output is desired (i.e., OPTION 6 is chosen), the order of 
input of the entries need not follow the above example. If the r-matrix is 
to be output, however, the program user must be careful of the order of his 
input data proportion sets, since their position in the r-matrix is determined 
entirely by their order in the input data. 

If an end-of-file mark is encountered before a complete matrix of [(K) (K-Dl/B 
Type III data sets is read, K is arbitrarily truncated to one more than the 
smallest value for which a square data matrix exists in the data. This will 
permit processing of all the proportional data, but will result in the use of 
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s 

indefinite operands if an attempt is made to output the r--matrix. Therefore, 
if proportional data (Type III) is used as input and the r-tnatrix is required 
as output, the user must be sure he includes exactly [(K) (K''l)]/2 sets of 
proportions, and incorporates them into his data in the correct order. If the 
r-matrix output is not required, any overestimate of NVAR-K will suffice. 

Omitting the READ INPUT DATA card results in the specified analyses being 
performed on whatever data (if any) are stored on file TAPEl. For the initial 
problem in a set of one or more problems, TAPEl is an empty file until the 
READ INPUT DATA card is vead. (A "problem" is defined as the analysis specified 
in the set of control cards between the PROBLEM card and the MATRIX card 
inclusive. More than one "problem" may be included in a single computer run.) 
After the initial problem, however, TAPEl contains whatever data was input 
for the preceding problem via the INPUT FORMAT card. Therefore, specifying 
a "problem" without including the READ INPUT DATA card, results in a second 
pass through the input data of the preceding problem. By varying NVAR, NO. 
OF CASES, and the INPUT FORMAT the program user can perform analyses of any 
subset of the initial problem^ s data, if the READ INPUT DATA card is omitted. 

I. The MATRIX Card (mandatory) 

Three options may be specified using the MATRIX card. The specification 
must be an integer number between 1 and 3, in the following format. 



flATRIX 1 For '^yP^ data, print the input proportions, tetrachoric 

li > > 4l5 1 T iUittitniqH»itf!fiiH:t and phi correlations, one variable pair at a time. Output 

the r-matrix in lower diagonal form. 

2 Output r'-matrlx in lower diagonal form only. (Suppress 
Individual variable pair results if Type III data are Involved.) 

3 Output r-matrix in square form. (Suppress individual variable 
pair results if Type III data are Involved.) 

The MATRIX card must be the last card in a problem set (excepting the 
FINISH card), since execution of the data processing and output options is 
begun after the MATRIX card is read. Placing the MATRIX card out of place 
will result in all cards following it being read as control cards for the 
succeeding problem or as data cards. Placing it before the READ INPUT DATA 
card will result in the analyses specified in the present problem being 
performed on the data of the preceding problem, or on an empty data file, as 
described above. 

If no matrix output is desired (i.e., OPTION 6 is elected) any MATRIX 
option (1,2 or 3) may be used to signal the start of data processing. 

J> llie FINISH Card (optional) 

This card has no specification field, and takes the form: 



3 <h I 7 tit till 1>l ? 3 W ' tsnh tiit»bin23?4b»17 i fe'3 i l ' it " ^>l»34y fi 
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When the FINISH card is encountered on the INPUT file, processing 
terminates, and control is returned to the computer operating system. Thus 
the FINISH card should appear only as the last card in a set of program 
control cards; it may not be used to separate problems. Failure to include a 
FINISH card results in printing an error diagnostic, program termination, and 
return of control to the computer operating system. In general, this will 
not affect the outcome of the user's programs other than TETREST, so the 
FINISH card may be considered optional. Its use is recommended whenever any 
data follow the TETREST control cards on the INPUT file. 

K. Control Card Order 

The following are the only restrictions on control card order within a 
problem. (A problem is the execution of the program on a single set or matrix 
of data. Multiple problems may be included in a single run of the program.) 
Aside from these restrictions, control cards may be placed in any order* 

1. The READ INPUT DATA card, if present, must follow the INPUT MEDIUM, 
INPUT FORMAT and NO. OF CASES cards whenever any or all are used. 

2. Data cards, if used, must follow the READ INPUT DATA card* If the 
input medium is CARD, and if the program is to count observations 
(number of cases) , then an end-of-record card (7-8-9) must be inserted 
between the last data card and the next program control card. Failure 
to comply with this requirement will result in program control cards 
being read as data cards, and will usually terminate program execution 
without a diagnostic. 

On the other hand, if the number of cases is specified exactly, no 
end-of-record card should separate the data cards from the control 
cards. 

3. The MATRIX card must be the last control card in each problem. 

4. The FINISH card must be the last card in a run. 
L. Examples of Control Card Use 

1. Card input of Type I data; program counts number of cases; a lower 
diagonal tetrachoric correlation matrix will be output. 
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RUN NAME EXAMPLE 1 

PROBLEM NVAR=05 R=TET 

INPUT FORMAT (5A1) 

OPTIONS 3 

READ INPUT DATA 

iOlOl 

illOO 



• • • • • 

00111 

> (RECORD SEPARATOR) 

MATRIX 2 
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2. Card Input of Type I data; number of cases specified exactly by 
user; a square phi correlation matrix will be output. 

RUN NAME EXAMPLE 2 

PROBLEM NVARS05 d=pht 

JNPUT FORMAT (SAD 

NO. OF CASES 0005 

OPTIONS 3 

READ INPUT DATA 

10101 

11100 

10111 

01100 

10110 

MATRIX 3 



3. Tape input of data; two problems in succession (one on cases 1-100 
in the data file; the second on cases 101-200); both tetrachoric and 
phi matrices will be output in lower diagonal format. 



RUN NAME FIRST 100 CASES 

PROBLEM NVARsOS R=BOTH 

INPUT MEDIUM TAPE 

INPUT FORMAT FIXED 

INPUT FORMAT (SAD 

NO* OF CASES 0100 

OPTIONS 2 

READ INPUT DATA 

MATRIX 2 

RUN NAME SECOND lOO CASES 

PROBLEM NVAR=05 R=BOTH 

INPUT MEDIUM TAPE 

NO* OF CASES 0100 

OPTIONS 2 

READ INPUT DATA 

MATRIX 2 

FI SK 



4. The above three examples are merely suggestive of the flexibility of 
options available to the program user. Far too many different con- 
figurations of control cards are permissible to list them all here. 
The following example uses the minimum number of control cards, 
setting pro'gram parameters by default everywhere possible: 



PROBLEM NVAR=05 
INPUT FORMAT {5X»5A1) 
READ INPUT DATA 

INSERT DATA HERE 
> (RECORD SEPARATOR) 

MATRIX 2 
FINISH 
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This minimal configuration assumes Type I input data punched on cards, 
with the program counting number of cases. It will print the 
tetrachoric correlation matrix only, in lower diagonal form and 8F10.7 
format, on the file OUTPUT. 

Additional examples of program input, in conjunction with the output 
they generated, are shown in Appendix A* 

Program Specifications 

File Structure 

Eleven files are used within TETREST, in the following order: 

INPUT 

OUTPUT 

PUNCH 

DATA the user-supplied file containing the data to be 

analyzed whenever the INPUT MEDIUM is DISC or TAPE. 

TAPEl the file used within the program to store the user's 

data. When the READ INPUT DATA control card is read, 
TAPEl is rewound, and data read from file INPUT or DATA 
is written onto file TAPEl. 

TAPES the file on which the r-matrix is written when the user 

chooses OPTIONS 4 or 5. TAPE 8 is never rewound within 
the program; an end-of-file is written on it after each 
matrix is written on it. The user must rewind TAPES 
in order to have access to the data on it. 

TAPE9 the file on which the joint frequency matrix is written 

when the user specifies DATA-2 or 4 on the PROBLEM card. 
The comments above for TAPES apply to TAPE9 as well. 

TAPE6-0UTPUT TAPE6 is a logical unit, set equivalent to the file OUTPUT. 

TAPEIO-PUNCH TAPEIO is a logical unit, set equivalent to the file PUNCH. 

TAPE51«DATA TAPE51 is a logical unit, set equivalent to the file DATA. 

TAPE55-INPUT TAPE55 is a logical unit, set equivalent to the file INPUT. 

Subprograms 

TETREST Includes the following subroutines and functions as part of the 
program: 

TETREST The main program; executes the input and output options. 
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SETUP Reads the program control cards and the input data; 

determines the desired input/output options; Issues 
error diagnostics* 

FREQl Constructs the KxK joint frequency matrix when Type I 

data are used. 

CELLS2 Constructs a 2x2 contingency table for each pair of 

variables; selects the appropriate frequencies for 
the analysis; constructs the rmatrlx. 

APXTET3 Approximates tetrachoric r by means of iteration given 

the joint and marginal proportions. 

GAUSS8 (Function) Performs 8-polnt Gaussian quadrature. 

MTRXOUT Outputs the appropriate correlation matrices In square 

or lover diagonal form. 



The following system library subroutines are used within TETREST: 



DATE 


subrcutiae returns the current day, month and yiar. 


CDFNI 


cumulative normal distribution function inverse. 


EXP 


the exponential function. 


AfiS 


absolute value function. 


SQRT 


the square foot function. 



Additionally, TETREST makes use of the Control Data ENCODE and DECODE 
statements, as well as the IF (EOF) end-of-file check. 

Timing and Core Memory Requirements 

TETREST requires 46000 rcLal core memory (CM) cells to run when using the 
FORTRAN extented (FTN) compiler under KRONOS 2.1 on the CDC CYBER 74 computer. 

Central processor (CP) time for a TETREST run is dependent on the data type, 
the correlation coefficient chosen, and the number of variables. 

Type I data requires the most CP time; Type III requires the least. For 
Type I data, CP time varies directly with the NO. OF CASES and with NVAR, the 
number of variables. 01 tiai^ for the other data types is a direct function 
of NVAR only. 

No formula is available for forecasting CP time requirements. Actual 
data runs have required CP time as listed below; 



Case Data Type NVAR NO. OF CASES CP Seconds 

I 75 80 15 

2 I 60 200 15 

3 I 50 240 25 

4 I 20 200 7 
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The disparity in CP time between Case 3 and Cases 1 and 2 illustrates the 
difficulty of prescribing a rule of thumb. CP time also depends on 
idiosyncracies in the data which affect the average number of iterations 
required per correlation. 

Availability 

A complete FORTRAN listing of the program is included in Appendix B. 
Card decks are available from the authors. A copy of the source code will 
be provided on a tape submitted by the prospective user. 



Method of Calculation 
Tetrachoric Correlation Coefficient 

The tetrachoric correlation between two dichotomous variables involves 
solution of the following equation for r, given the values of L, h,k,x, and y. 



L(x,y,r)« 



exp 



■(h^+k^-2rhk) 
2(l-r^) 



dh dk 



(1) 



where r is the tetrachoric correlation 
x,y are standard normal deviates 

L is a joint function of the two dichotomous variables 
and their tetrachoric relationship (whose numerical 
value is observable as a joint proportion) 

The values of x and xjy h and k are obtained from calculations based on the 
marginal proportions of the two dichotomous variables, which can be obtained 
in turn from their 2x2 contingency table: 
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Variable 1 
incorrect correct 





(0) 




Correct (1) 

Variable 


b 


a 


i 

Incorrect (0) 


d 


c 




b+d 


a+c 


proportion 







24 



total proportion 



a+b 



c+d 
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The parameters of Equation 1 are estimated from the sample contingency table as 
follows : 

x«*(q^)*the standard normal deviate corresponding to q^(*l-p^) 
y«*(q^)"the standard normal deviate corresponding to qj(«=l-pj) 

L=»d/N«the joint proportion of persons responding positively to both 
variables (items) 

In general, approaches to solving Equation 1 have involved the use of 
an infinite series • McNemar (1955) gives one such solution, expressed here 
in notation slightly different from his: 



hk 

where h,k 



N ' ^i^j ^ rfxy r' 



2! 



. ^ 2 3 
(x-'-l) (y -l)_rl 



(2) 



3t 



are the ordinates of the standard normal curve at 
X and y, respectively 



McNemar •s formulation includes the restriction that q^ and q^ (the marginal 

proportions) must each be less than or equal to 1/2 • Obviously, that 
restriction implies that t may be the joint frequency value occurring in cell 
ay by Oy or d of the 2x2 table, whichever cell corresponds to the two marginal 
values meeting the restriction. Use of an inappropriate cell to obtain tTie 
joint frequency t will introduce an error in the proportion t/N and hence will 
in general result in an incorrect numerical solution for v. TETREST, in 
dealing with Type I data, constructs the 2x2 table for each pair of variables 
and selects the appropriate joint frequency for correct solution of Equation 2. 

Kirk (1971, 1973) devised an approximation procedure based on Gaussian 
quadrature and use of a Newton-Raphson iteration procedure to evaluate the 
integral and approximate r». The first two terms of Equation 2 are used to 
obtain an initial estimate based on the empirical values of t/N, q^ and q^ • 

In order to employ Gaussian quadrature and restrict the range of the 
estimate of r». Kirk gives the following transformation, based on Equation 1. 



f(r)= 



27r 



/1-u r 



exp 



2 2 
■(x -2xyur-Hy ) 



2(l-u^r^) 



du 



(3) 



27r 



i=0 



w^g(u^) 



(4) 



where g(u) 



/1-u r 



exp 



2 2 
-(x -2xyur4y ) 



2(l-u^r^ 
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and the are the roots of the Legendre polynomials; the are the associated 
weights for an (n+l)^poinc quadrature. 

The starting value for solution of Equation 4 is obtained by direct 
solution of the following equation for r^: 

^^"^^ . r^+xy rj (5) 



hk 



2i 



where P«t/N, taken from the appropriate cell of the 2x2 table. A limiting 
value of +«80 is used whenever the absolute value of the starting estimate 
exceeds 1.00. Successive values are calculated using Newton-Raphson iteration: 



r » T -• 



(6) 



f (r^) 



where in-P-qj^q^ 
and 



f ' (r) - exp 



2 2 
r- (x -2xyr+y 1 

L 2(l-r^) J 



(7) 



Iteration continues through as many as 100 cycles or until 

II /a\ 
'^i+l"'^ii <.0001 

Convergence failures can occur. When convergence to the criterion of .0001 
does not occur, one of two alternatives takes place: 

1. If convergence to within .001 has occurred, an estimate of r is 
returned, but an informative diagnostic is printed, telling the user 
the indices of the two variables involved and the degree of convergence 
attained. 

2. If convergence to within .001 has not occurred, a value of +99.00 
is returned for 

Generally, failures of convergence are due to extremely low values or 

the joint proportion, P, or to the joint proportion P being very close to 

either of the marginal proportions, ox q.* One circumstance in which this 

^ Q 

is a problem is the case in which one or more cells in the 2x2 table contains 
a zero. 
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When three zero cells occur, as In the following 



a 


0 


0 


0 



the tetrachoric correlation is automatically set to zero, since there is no 

covariance. Similarly, when there is zero variance in only one 

variable (p^^-LOO or Pj=»1.00) there is zero covariance between variables, 

and the tetrachoric correlation ?.s zero. 

Whex^ just two zero cells occur, and these occupy either diagonal of 
the 2x2 table, the program returns a value of +1.00, the sign being determined 
by the sign of the covariance between the two variables. For example: 



j 



0 


b 


c 


0 



here the covariance equals j > correlation coefficient is returned 

as -1.0. 



In the case of a single zero cell, such as 



1^0 



0 


b 


c 

b- — — 


d 



there is no rational basis for setting the tetrachoric r to any fixed value 
as was done above. There is non-zero covariance, as well as variance in 
both variables. This class of problems frequently results in failure of the 
iteration procedure to converge on a value of the tetrachoric r. Non-convergence 
is especially a problem with a single zero cell concurrent with an extreme 
proportion in both variables. As indicated earlier, a failure of convergence 
is signalled by a value of +99.0 being returned in place of the tetrachoric r. 

Kirk's algorithm for approximating the value of the tetiachorlc correlation 
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does not directly calculate the correct sign of the relationship. For Type I 
and Type II data, the sign of the tetrachoric correlation is determined by 
the sign of the covariance, which is completely determined by the determinant 
ad-bc of the fourfold table of any variable pair. Thus, when ad-^ba is positive, 
a positive value for the tetrachoric is output; a negative value is output 
when ad-be is negative. 

For Type III data the sign of the tetrachoric r is corrected by the 
sign-change character, SIGN. If SIGN was blank on input, the value calculated 
by TETREST is returned uncorrected. If SIGN was any non-blank character, 
the calculated value of the tetrachoric v is multipled by -1 to correct the 
sign of the relationship. 

The Phi Coefficient 

For Type I and II data, the value of the phi coefficient is calculated 
from the frequencies of the fourfold table, using the following formula: 

PUI . ad-bc (9) 

/(a+b) (a+c) (b+c) (b+d) 

In the case of Type III data, phi is calculated from the input proportions 
and SIGN: 

PHI - - "i'^j (10) 

/(qi)(i-qi)(qj)(i-qj) 

where SIGN ■ +1 or -1 according to the input convention described above. 
Accuracy 

The method used for estimating the tetrachoric correlation coefficients 
in TETREST is a numerical approximation technique. Kirk (1973) reported 
results generally accurate to three decimal places, with few exceptions, using 
his algorithm and double-precision arithmetic. TETREST involves a slight 
modification of Kirk's method, and does not use double-precision. Nevertheless, 
on a CDC CYBER 74 computer, TETREST results are still usually accurate to 
three decimal places. Users may evaluate for themselves the accuracy of 
TETREST by using, as Type III input data, joint and marginal proportions 
corresponding to known true values of the tetrachoric correlation. National 
Bureau of Standards (1959) tables are one source of such data. 

The phi coefficient is calculated exactly, and is accurate to 15 significant 
digits on Control Data machines. 
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Appendix A 
Examples of Input and Output 



PROBLEM NVAR=09 DATA=a 2X2= NFMT= key-1 R-rOTH 

RUN NAME DEMO-TYPE I DATA INPUT ON FILE INPUT 

INPUT FORMAT (20X?9A1////) 
OUTPUT FORMAT (1X»20F6.2) 
READ INPUT DATA 

END-OF-PILE ENCOUNTERED WHILE READING DATA 
NO, OF CASES HAS BEEN CHANGED TO 66 
MATRIX 2 



DEMO—TYPE I DATA INPUT ON FILE JNPUT 
PRI».TED BELOW IS THE JOINT FREQUENCY MATRIX* 

WITH JOINT FREQUENCY CORRECT»++» IN THE UPPER DIAGONAL* 

JOINT FREQUENCY INCORRECT »— » IN THE LOWER DIAGONAL 
AND MARGINAL FREQUENCY* + ,IN THE DIAGONAL 
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uisted below is The matrix of tetrachoric correlations 



VARIABLE 1 

VAN TABLE 2 
.69 1.00 

VARIABLf 3 

.73 t^^S 1.00 

VARIABLE ^ 

-.05 -.05 1.00 

VARIABLf 5 

.29 .28 .**9 .29 1.00 

VARIABLE 6 

•56 .59 .24 .28 .63 1.00 

VARIABLE 7 

.50 .69 .19 .53 -.07 .55 1.00 

VARIABLE 8 

.75 .77 .50 .45 .06 #62 .77 l.OO 

VARIABLE 9 

.76 .76 .67 ,44 ,73 .77 .81 .76 1.00 

LTfiTED BELOW IS THE MATRIX OF PHI CORRELATIONS 



VARIABLE 1 
1.00 

VARIABLE 2 
.42 1.00 

VARIABLE 3 . 
.46 .29 1.00 

VARIABLE 4 , 
-.02 .28 -.03 1.00 

VARIABLE 5 

.16 .18 .32 .18 1.00 

VARIABLE 6 

.35 .**0 .15 .18 .43 1.00 

VARIABLE 7 

.27 .'♦7 .12 .35 -.0** .35 1.00 

VARIABLE 8 

•52 •52 •32 ^27 .04 .41 •^S 1.00 

VARIABLE 9 

.52 .53 .**6 .27 .51 .55 .51 .55 1.00 
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LISTED BELOW IS THE MATRIX OF TETRACHORIC CORRELATIONS 



VARIABLE 1 
l.OOOOUOO 

VARIABLE 2 
.<^9q9974 l.OOOOOOO 

VA,<IABLE 3 
«BOO00o0 .8500160 I.QOOOOOO 

VARIABLE 4 
«A486960 .8999555 .9000036 l.OOOOOOO 

VAwIABLE 5 
.q50l760 .8985307 .9500787 .9601756 l.OOOOOOO 

VARIABLE 6 

,9601208 .9910074 .9910019 .9910300 .9910300 1.0000000 
VAkIABLE 7 

-.0500211 -.1994941 -.7499432 -.7500099 -.9000033 -.9017240 l.OOOOOOO 
VARIABLE 8 

-.9500391 -.9600849 -.9597555 -.9910023 -.9571481 -.9701562 1.0000000 
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Appendix B 

A Complete FORTRAN Listing of PROGRAM TETREST 



PROGRAM TETREST ( INPUT » OUTPUT » PUNCH » TAPE 1 » TAPES t TAPE9 » DAT A » TAPE6=0U I 

♦TPUT f TAPEl 0=PUNCH » TAPE5l=DATA » TApE55=INPUT ) 2 

C0MM0N/RFREQ/R(75»75) ^ 

COMMON/FRMAT/IFMTiaO) »0FMT(8) ^ 

C0MM0N/I0PT/0PT(8) = 

COMMON/MISC/ ICOL ( 75 ) » NAME ( 8 ) » N I » nS » NF » I R »L lU » JSPEC 6 

DIMENSION DUM<75) ' 

INTEGER OFMT » 

IRUN=0 9 

1 CALL SETUP 1? 

CALL OATE(NOw) 1* 

PRINT 1019 » NAME* NOW 12 

ZL=NS 13 
REWIND 1 

IF(0PT(1) .EQ.5.)G0 TO 7 15 

IF{0PT(1).GE.3.)60 TO 31 16 

CALL FREQ1{NS»NI»IR) 1' 

IF(0PT(1) .EQ.2.)U»61 1® 

31 READ(1,IFMT) (<R(X»J)»J=1»NI),I=1,NI) 19 

IF(0PT(1) .EQ.«l.)'»»61 20 

H PRINT 1003 21 

DO 5 1=1 »NI 22 

5 IC0L(I)=I 23 
PRINT 100«*» (ICOL(I) »I=1'NI) 2H 
PRINT 1005 25 
DO 6I=1»NI 26 
PRINT 1006»I 27 
ENCODE(l25»1007»DUMtl))<R(I»J) »J=1»NI) 28 
DECODE<125»1008»DUM<1))(ICOL(J)»J=1»NI) 29 

DO 51 J=1»NI 5? 

IF(IC0L(J).EQ.'*H )IC0L(J)=4H 0 31 

51 IC0L(J)=IC0L(J) .A.77777777007777777777B.O.550000000000B 32 

IF(0PT(1) .NE.U.)GO TO 6 33 

WRITE(9»1008) (ICOL(J) »J=1»NI) 3H 

6 PRINT 1009» (IC0L(J)»J=1'NI) 35 
ENDFILE 9 3° 

61 CALL CELLS2 H 

GO TO 10 38 

7 1^=0 S 
IF(0PT(2) .GE.2.)eO TO 71 

PRINT 1012 



71 



72 



PRINT 1013 



ASI6N=1. 



HZ 



DO 8 I=2»NI 
R(I»I)=1. 
K=I-1 

DO 8 J=1»K 

R(I»J)=R(J»I)=99.0 ^, !J| 

READ (1»IFMT)PC»PI»PJ»ISIGN '♦j 



«»5 



50 



iFdSlGN.NE.lOH )ASI6N=-1. 5l 



52 
53 

55 



ERLC 



IF(E0K(1) )99U»72»99U _- 
CALL APXTET3(PC»PI»PJ»0UT»EPS) |J 
IF(ASI6N.EQ»0.)ASIGN=1. 
R(J»I)=OUT*ASIGn 
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PHI = (PC-{PI*PJ) )/i,QRT(Pl*PJ*(l.-pI)*(l.-pj) ) 


56 




PHI=PHI*ASIUN 


57 




IF(0PT(2) .GE.2.)G0 TO 73 


58 




0UT=R(J»I) 


59 




N=N+1 


60 




PRINT 10m»N»I»J»PC»PI »PJ»PHI»OUT»EPS 


61 


73 


CONTINUE 


62 


8 


RU»J)=PHI 


53 




IF(0PT{2) .GE.2.)60 TO lO 


64 




PRINT 1010 


6S 


CCCCC 


SUPPRESS R-MATRIX 


66 


10 


IF{0PT(3) .EQ.6.)G0 TO 1 


67 




ISPEC=1 


68 




IFIJSPEC,EQ.3HB0T) ISPEC=2 


69 




IF(JSPEC,EQ.3HPh1)ISPEC=.^ 


70 


CCCCC 


BRANCH TO R-MATRIX OUTPUT 6R0UP 


71 




IF(0PT(3) .LE.i^OlOO^llO 


72 


100 


IF{0PT(3) .EQ.DGO TO 105 


73 




LU=6 


7** 




IF(0PT(2) .EQ.3.)10l»10i> 


75 


101 


ASSIGN 102 TO NSTATE 


76 




GO TO(20l»201»202) »ISPEC 


77 




WRITE<6»9001) 


78 




GO TO 1 


79 


102 


IF(0PT(3) .EQ.3. )(50 TO 1 


80 




IF (OPT (3) .EQ.2.)105»110 


81 


CCCCC 


print»lower diagonal 


82 


103 


ASSIGN 104 TO NSTATE 


83 




GO To (203»203»20'*) »ISPEC 


8*^ 




mRITE(6»9001) 


85 




GO TO 1 


86 


lOH 


IF(0PT(3) .EQ.3.)G0 TO 1 


87 




IF(0PT(3) .EQ.'*.)GO TO 110 


88 


CCCCC 


PUNCH SQUARE MATRICES 


89 


105 


ASSIGN 108 TO NSTATE 


90 




Luno 


91 




IF(0PT(2) .E'3.3.)106»107 


92 


106 


GO TO(20l»20l»202) 'ISPEC 


93 




WRITe(6»9001) 


9<» 


107 


60 TO (203»2o3»20U)»ISPEC 


95 




«|RITE(6»9001) 


96 


108 


GO TO 1 


97 


CCCCC 


WRITE R-MATRICES 


98 


110 


IF(0PT(2) .EQ.3.)lll»113 


99 


CCCCC 


SQUARE MATRICES 


100 


111 


IF(ISPEC.EQ.3)G0 TO 112 


101 


CCCCC 


SQUARE TET 


102 




CALL MTRX0UT{8»1) 


103 




IF(ISPEC.EQ.2)112»1 


10<I 


CCCCC 


SQUARE PHI 


105 


112 


CALL MTRX0UT(8»2) 


106 




GO TO 1 


107 


CCCCC 


LOWER DIAGONAL MATRICES 


106 


113 


IF(ISPEC.EQ.3)G0 TO 11*+ 


109 




CALL MTRX0UT(8»t*) 


110 
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IF(ISPEC.EQ.2)11(+»1 m 

im CALL MTRX0UT{8»3) ll2 

GO To 1 113 

CCCCC SQUARE. TET H** 

201 CALL MTRX0UT(LU»1) JlS 
IF(ISPEC.EQ.2)202»2021 Jl^ 

CCCCC SQUARE PHI H' 

202 CALL MTRX0UT(LU»2) 

GO TO 2021 119 

203 CALL M7RX0UT(LU»U) 120 
CCCCC LOWER TtT 121 

IF(ISPEC.EQ.2)20U»2021 122 

CCCCC LOWER PHI 123 

20'* CALL MTRX0UT(LU»3) 12^ 

2021 GO To NSTATE» (102»iOU»l08) 125 

GO TO 1 126 

1001 FORMAT(8A10) 127 

1002 FORMAT(lX»aA10) 128 

1003 F0RMAT{1X»*PRINTED BELOW IS THE JOINT FREQUENCY MaTRIX»WITH JOINT 129 

2 FREQUENCY CORRECT»++» IN THE UPPeR DIAGONAL* */51X»*J0INT FREQUENCY 130 

3 INCORRECT » — f IN THE LOWER Dl AC-ONAL* »/5lX » *AND MARGINAL FREQUEN 131 
UCY» + »IN THE DIAGONAL*/) 132 

lOO** F0RMAr(7x»26I5) 13^ 

1005 FORMAT(//) 13^ 

1006 F0RMAT{lX»i3»X»lH*) 135 

1007 FORMAT(25F5.0) 136 
iOOB FORMAT ( 25 (AU» IX ) ) 13; 

1009 F0RMAT(7X»25A5) 138 

1010 FORMAT(lHl) ^ ^„ 139 

1011 F0RMAT(1X»*PRINTED BELOW ARE THE 2X2 CONTINGENCY TABLES FOR EVE UO 
2RY item pair*//) , l**! 

1012 FORMAT ( IX » ♦THE FOLLOWING TETRAcHORIC CORRELATIONS WERE CALCULATE U2 

2D FROM THE InPUT PROPORTIONS LISTED.* »/lX»*CAUTlON ANY MARGINAL 1U3 

3FREQUENCIES GREATER THAN .50 WERE CORRECTED INTERNALLY* WHICH CAUSE 

US AN INCORRECT RESULT* »/lX»*NOTE THE SIGN OF THE CORRELATION 1U5 

5C0EFFICIENT MAY BE REVERSED IF ThE JOINT PROPORTION WAS TAKEN FROM U6 
6 THE MINOR DIAGONAL**/) 

1013 FORMAT(1X»7HPROBLEM»3X'9HVARIABLES»5X»30HJOINT AND MARGINAL PROPO U8 
2RTI0nS»16X»3mPHI»7X»11HTETRACH0RIC»5X»7HEPSIL0N/) 1^9 

lOlU FORMAT(lX»lb»5x»l«*»*»*»I'*»X»3(Fi0.8»5X) »3(F10.5»5X) ) 150 

1015 F0RMAT(/1X»*VARIABLE*»I5) 151 

1017 FORMAT(8F10.7) }|; 

1018 F0RMAT(10F8.5) 

1019 FORMAT(1H1»8A10»UOX»A10/) JSH 
9001 fORMAT(/10(1H*) »*PROGRAM ERROR ISPEC IS OUT OF RANGE INTETReST*) 155 

END 



SUBROUTINE MTRXOUT (LU» iSPEC) JS^ 

CCCCC OUTPUTS THE UPPER OR LOWER DIAGONAL OF A SQUARE MATRIX AS A SQ 158 

CCCCC OR LOWER DIAGONAL MATRIX }=J 

C0MM0N/RFREQ/R(7b»75) J^" 

COMMON/FRMAT/IFMT(80) »t.-MT(8) 
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COMMON/MISC/DUM(75) »NAME (8 ) »NI »NS 'NF » JR » LIU* JSPEC 162 

1F(LU.E«.6)99»100 163 

99 GO T0(105.106»106»105) »ISPEC le^ 

100 60 T0(l02,101»l0i+»103) f ISPEC 165 
CCCCC.»» ••LOWER DIAGONAL TO SQUARE 166 

101 DO H I=1»NI 167 
IF(LU*EQ.6)1»2 168 

1 WRITE<6'1015)I 169 

2 iFd.EQ.NDGO TO b 170 
K=I+1 171 
DO 3 J=K»NI 172 

3 DUM(J)=R(K»1) 173 
H WRITe(LU»OFMT) ( (R(I»J) »J=1»NI) , (dUM(L) »L=K»NI) ) 

5 WRITe(LU»OFMT) (R(l»J} » J=1»NI) 175 
GO TO 20 176 

CCCCC UPPER DIAGONAL TO SQUARE 177 

102 IF(LU.EQ.6)6»7 178 

6 1=1 179 
WRITE(6»1015)I 180 

7 WRITE(LU»OFMT) (R(1»J) »J=1»NI) I8l 
OO 9 I=2»NI 182 
IF(LU.EQ.6)71»72 183 

71 WRITE(6»1015)I 18** 

72 IFd.EQ.NDGO TO 10 185 
DO 8 J=1»NI 186 

8 DUM(j)=R(J»i) 187 
L=I+1 188 

9 WRITE(LU»OFMT) ( (DUM(J) »J=1»I) » (R(I»K) ,K=L»NI) ) 189 

10 WRITE(LU»0FMT) (R(I»NI) »I=1»NI) 190 
GO TO 20 191 

CCCCC UPPER DIAGONAL TO LOWER 192 

103 DO 13 I=1»NI 193 
DO 11 J=1»I 19'* 

11 DUM(J)=R(J»I ) 195 
IF(LU.EQ.6)12»13 196 

12 WRITE(6»1015)I 197 

13 WRITE«LU»0FMT) (DUM(J) »J=1»I) 198 
GO TO 20 199 

1015 F0RMAT(/1X»*VAR1ABLE*»I5) ?00 

CCCCC. 'LOWER DIAGONAL TO LOWER 201 

lOH DO 15 1=1 »NI 202 

K=I 203 

IF(LU.EQ.6)1'*»15 20'» 

11 WRITE(6»1015)I 205 

15 WRITE(LU»OFMT) (R(i»J) » J=1»K) 206 

GO TO 20 207 

105 WRITE(6»1051) 208 
GO TO 100 209 

106 wRITE(6'l06i) 210 
GO T0(102,101»10'*»103) »ISPEC 211 

20 END FILE LU 212 

RETURN 213 
1051 F0RMAT(*1LISTED BELOW IS THE MATRIX qF TETRACHORIC CORRELATIONS* 211 

+///) 215 

1061 F0RMAT(*1LISTED aELOW IS THE MATRIX OF PHI CORRELATIONS*///) 216 
END 
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SUBROUTINE SETUP flj 

C0MM0N/RFREQ/R(75»75) • 219 

COMMON/FRMAT/IFMT(80) »0FMT(8) 220 

C0MM0N/I0PT/0PT(8) "1 

COMMON/MI SC/I COL ( 75 ) »NAME (8)»nI»NS»NF»IR»LIU» JSPEC 222 

DIMENSION IN(8)»DUM(75) 223 

NS=9999SNF=1SLIU=55$NI=0 22«J 



PRINT 2007 



225 



CALL DATE (NOW) 226 



IVER=9H 1 FEB 75 



227 



PRINT 1000»NOW»IVER |2f 
DO 2 1=1 »8 ItH 



NAME<I)=10H 
0PT(I)=I 
0FMT(I)=10H 



230 

231 
232 



OFMT{l)=10H(aFl0.7) ?33 



11 DO 12 I=l»8 



23** 
235 



12 IN(I)=10H 

READ(55»2001)IN f^; 

IF(E0F(55))991»13»991 |37 
)3 PRINT 2002»IN 

iFdNd) .EQ.IOHPROBLEM )60 TO 20 23^ 

iFdNd) .EQ.IOHRUN NAME )60 TO 30 2JjO 

IrdNd) .EQ»10HINPUT MEDI)60 TO UO 2i»l 

IFdNd) .EQ.IOHINPUT F0RM)60 TO 50 2<t2 

IFdNd) .EQ.IOHNO. OF CAS)60 TO 60 2^3 
IFdNd) .EQ.IOHOUTPUT F0R)60 TO 70 

IFdNd) .EQ.IOHOPTIONS )50 TO bO 2Hb 

IFdNd) .EQ.iOHMATRIX )eO TO 9O 2Hb 

IFdNd) .EQ.lOHREAD INPUT)60 TO lOO i^' 

IFdNd) .EQ.IOHFINISH )G0 TO 110 |J° 

PRINT 992 iZf, 
ST0P2 

20 DECODE(60»2003dN(3))Nl»OPT<l)»OpT(U) ,NF»IR» JSPEC 251 

IF (JSPEC .EQ.3H )JSPEC=3HTET 252 

lF(NI.LE.1.0R.NI.6T.75)eO TO 993 1=3 

lF(OPTd).LT.0..OR.OPT(l).6T.5.)GO TO 993 25H 

iFCOPTd) ,EQ.O.)OPT(1)=1. 255 
iFdR.EQ.SH )IR=1H1 
IF(NF.EQ.0)NF=1 

IF(NF.LT.0. OR. NF.tiT. 10)00 TO 993 

IF(OPT(U).LT.0..OR.OPT('*).6T.9.)gO TO 993 259 

60 TO 11 m 

30 DO 31 I=l»7 

31 NAME(I)=IN(I+1) 2°3 
nAME(8)=10H 26U 



60 TO 11 



265 



j»0 IF(IN{2).EQ.10HUM CARD )UliU2 " 
ifl LIU=55 * 60 TO 11 2ft7 
<f2 iF(lN(2).EQ.10HUM TAPE )U3»UU 



lis 



-34- 



U3 LIU=51 S GO TO 11 268 

IF(IN(2) .EQ.IOHUM DISC )G0 TO 269 

PRINT 9009 270 

GO TO 11 271 

bO DO 51 I=l»80 272 

51 IFMT(I)=10H 273 
DEC0dE(70»200**»IN(2) ) ( IFMT ( I ) » 1=1 »7) 27*+ 
IF(NF.EQ.I) H»52 275 

52 nF=7+((NF-1)*8) 276 
READ(55»200l) ( lFt.;T( I ) » I=8»NF) 277 
PRINT 2002» (IFMT(l) »I=8'NF) 278 
GO TO 11 279 

bO DEC0DE(l0»2005»Ih(2) )NS 280 

IF(NS.LT.1.0R,NS.GT.9999)G0 TO 99** 281 

GO TO 11 282 

70 DECOdE(70»200**»IN(2) ) (OFMTCI) »I=l»7) 283 

60 TO 11 28** 

80 DEC0dE(10»2006»IN(2) )0PT(3) 285 

IF(0PT(3) .LE.0..0R.,0PT<3) .GT.6.)G0 TO 995 286 

60 TO 11 287 

90 DECOdE(10»200&»IN(2) )0PT(2) 288 

IF(OPT(2).LE.O..OR.OPT(2) .6T.3.)60 TO 996 289 

RETURN 290 

100 REWIND 1 291 
IF(0PT(1).EQ.5.)101»103 292 

101 DO 102 I=2»NI 293 
N1=I-1 294 
DO 102 L=1»N1 295 
READ(LIU»IFMT) (DUM(J) »J=1»**) 296 
IF(E0F(LIU) )1021»102,1021 297 

102 WRiTEd'iFMT) (DUM(J) .J=l'**) 296 
GO TO 1022 299 

1021 NI=I 300 

1022 END FILE 1 l°t 
60 TO 11 fO* 

103 IF(OPT(1).GE.3.)10«*»106 ^03 
10<l DO 105 I=1»NI 

READ(LIU»IFMT) (DUM(J) f J=1»NI) 30b 

IF(E0F(LIU) )997»105»997 306 

105 WRITEd'IFMT) 'C(UM(J) »J=1'NI) ^0' 
END FILE 1 308 

60 TO 11 5?0 

106 DO 107 I=1»NS 

READ(LIU»IFMT) (IC0L(J)»J=1»NI) 311 

IF(E0F{LIU))1061»107»1061 ^1^ 

1061 NS=I-l flj 
PRINT 9008»NS 

GO TO 108 l\l 

107 wRITEd'iFMT) (ICOL(J) »J=1»NI) 31b 

108 END FILE 1 l^Z 
60 TO 11 ^18 

110 STOP tiZ 

991 CONTINUE :J|, 
PRINT 9001 

PRINT 9999 ^'^ 
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993 

99U 

995 

996 

997 

9001 
9003 

9004 
9005 
9006 
9008 
9009 
992 

9997 

9999 
1000 

2001 
2002 
2003 
2004 
2005 
2006 
2007 



9003 
9999 

9004 

<i99: 



9999 

9006 
9999 



rJOPl 
PRINT 
PRINT 
ST0P3 
PRINT 
PRINT 

ST0P4 
PRINT 9005 
PRINT 
ST0P5 
PRINT 
PRINT 
STOP 6 
PRINT 9997 

-•PRINT 9099 - - - ^ -- . _ „ - . 

STOP? 

F0RMAT(/10(1H*) .*ERROR— END/OF-FILE ENCOUNTERED WHILE READING TE 

2REST CONTROL CARDS*) . ^ „ .urr^ 

FORMAT ( /lO ( IH* )» TERROR— ONE OR MORE OF YOUR PROBLEM CARD PARAMET 
2ERS IS OUT OF RANGE*»/18X»*CHeCK THE vALuES OF NVAR»NF»DATA AND 
32X2.*) 

FORMAT ( /lO ( IH* )» ♦ERROR THE NO. OF CASES IS TOO LARGE OR NON-NUME 

2RIC VALUE WAS READ*) 

F0RMAT(/10tlH*) »*ERROR THE OPJION CARD VALUE IS OUT OF RANGE OR 

2 A NON-NUMERIC VALUE WAS READ*) 

FORMAT(/10(1H*)»*ERROR THE MaTRIx CARD VALUE IS OUT OF RANGE 0 

2H A NON-NUMERIC VALUE WAS READ*) , ^ , 

FORMAT{/10X»*END-OF-FILE ENCOUNTERED WHILE READING DATA*»/ 

2lOX»*NOvOF CASES HAS BEEN CHANGED TO *»I8) 

FORMAT (/10X»*COMMENT CHECK YOuR INPUT MEDIUM CARD FOR ERROR. IN 

2PUT MEDIUM IS ASSUMED TO BE CARDS**) 
FORMaT(/10(1H*) »*£RROR THE PROGRAM DID NOT READ A COMPLETE SET 

20f control cards. *»/18x»uf control cards are complete » check for m 
3:sspelled definition fields*) 

format ( 3 1 ( ih* ) »*error end-of-flle encountered while reading joi 

2nt frequency matrix*) 

format (/10x.*fijrther diagnostics would be misleading.*) 

FORMAT (100X»Alo/100X»*PROGRAM TeTREST*»/100X»A10»2X» *REVISION*» 
+//////////) 
F0RMAT(8A10) 
FORMAT{20X»8A10) 

FORMAT ( I2» 6X»F1»9X»F1»9X»I2» ax »Al»9X» A3) 

FORMAT(5X»A5»6A10) 

F0RMAT(5X,»I4»1X) 
F0RMAT(5X»Fi.4X) 
FORMAT (IHl) 
STOP 
END 



323 
324 
325 
326 
327 
328 
329 
330 

331 
3?? 
33- 
334 
335 
336 
-337 
338 
339 
340 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 

351 
352 
353 
354 
355 
356 
357 
358 
359 
360 

361 
362 
363 
364 
365 
366 
367 
368 
369 
370 



SUBROUTINE FREQl (NS»NV» IR) 
COMMori/RFREQ/IFR ( 75» 75) 
C0MM0N/FRMAT/IFMT(80) 



371 
372 
373 
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DIMEnSION IX(75) 374 

REAL IFR 375 

Nl=NV 376 

D0102 I=1»NI 377 

D0102 J=l»Nl 378 

IFR(I»J)=0. 379 

102 CONTINUE 380 
00 107 L=1»NS 381 

103 READ(1»IFMT) (IX(I)»I=1»NI) 382 
DO 107 I=lfNI 383 
DO 107 J=I»NI 38'* 

lF(IX(I).E(i.IX(J))lO'*»107 385 

104 iFdXd) .EQ.IR)lOb»i06 386 

105 IFR(I»J)=IFR(I»J)+1. 387 

- SC'-TG 107- - - - -. - . ™ , . 3458 - . 

106 IF(I.EQ.J)60 TO 107 389 
lFR(;j»I) = IFR(J»I)+l. 390 

107 CONTINUE 391 
RETURN 392 
END 393 



SUBROUTINE CtLLS2 394 

C0MM0N/RFreQ/R(75»75) 395 

COMMON/ I OPT/OPT (3) 396 

C0MM0N/MISC/IC0L(75) »NAME(8) »NI »n2»NF » IR »LIU » JSPEC 397 

Nl=NI-l 398 

IF(0PT(4).NE.9.)G0 TO 200 399 

PRINT 3 '♦OO 

3 FORMAT (*1PRINTED BELOW ARE THE 2X2 CONTINGENCY TABLES FOR EACH VA 401 
2RIABLE PAIR*»/1X»*99.00 IS PRINTED WHEN THE PROGRAM CANNOT COnVERG 402 

3E ON A SOLUTION**//) '♦03 

200 DO 227 1=1. Nl 

N2=I+1 '♦05 

DO 227 J=N2»NI 1*0^ 
A=P(I»J) 

U=R(J.I) '♦0^ 

E=R(I»I) '♦09 

G=R(J»J) '♦1? 



B=E-A 



F=C+D 



411 



C=G-A '♦1? 



413 



H=B+D '♦l^ 
ZL=A+B+C+D '♦^f 
DET=(A*D)-(B*C) 
IF(OET)2002»2001»2002 

2001 PHI=0. 

R(I»J/=0. l'-^ 
GO TO 2003 ^20 

2002 PHI=DET/SQRT(E*F*G*H) ^21 

2003 IF ( JSPEC. EQ.3HPHI)G0 TO 226 

lFiF.6E.E)201»202 ^f^ 
201 ANUMrE '♦2'* 
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IBITS2 $ GO To 203 "^SS 

202 ANUMsF '♦26 
iBITsO «*27 

203 lFlH,GE.G)20'+»205 '♦28 

204 BNUMsG '♦29 
IBIT=IBIT+1 $ 60 TO 206 '♦30 

205 BNUM=H '♦31 

206 IF(IbIT.EQ.O)207j208 '♦32 

207 CELL=D $ 60 TO 215 '♦S^ 

208 IF(IBIT.EQ.I)209»210 '♦34 

209 CELLsC S GO TO 215 '♦35 

210 IF(IBIT.EQ.2)211»212 '♦36 

211 CELLsB S GO TO 215 '♦37 

212 lF(lBlTtE0.3)214»2l3 '♦38 

213 . STOP200.0JL . .... „ _ . '♦.39 . 

214 CELLsA '♦40 

215 P=CELL/2L '♦41 
Ql=ANUM/ZL '♦'♦2 
QJsBNUM/ZL {^'♦3 

216 IF(DET)220»226»220 '♦'♦*♦ 
220 CALL APXTET3(P»QI»QJ»R1'EPS) '♦45 

IF ( DET ) 222 » 223*223 '♦'♦6 

222 X=-l. $ GO TO 225 '♦'♦" 

223 X=l. $ GO TO 225 '^'♦g 

225 R(I»J)=X*ABS(R1) '♦'♦f 

temp=abs(ri) '♦so 

lF(T£MP.EQt99t)60 TO 226 '♦Si 

IF(EPS.LE..0001)GO TO 226 '♦52 

PRINT 2»I»J»EPS '♦53 

226 R(J'I)=PHI '♦f^ 
IF(0PT(4).NE.9.)G0 TO 227 '♦SS 
PRINT 1»J»B'A»E»P»R(I»J) 'EPS»I»D»C»F»(5l»G»H»ZL»Qd»PHl '♦56 

227 CONTINUE '♦57 
DO 228 1=1 »NI ^SS 

228 R(I'I)=1. f59 

1 FORMAT (17X»6H- VAR»I3»3H ♦»8X»6HT0TaLS»11X»11HPR0P0RTI0NS/ 460 
212X»1H+»3F10.0»5X»6HP(I»J) »F9.8»9X»17HTETRACH0RIC R =»Fl3.7»3X» 461 
39HEPSIL0N =»F10.5/lX»3HVAR»I3»5X»lH-»3F10t0»5X»6HQ( I )»F9.8/ 462 
47X»6HTOTALS»3F10.0»5X»6HQ( J ) »F9.8»20X»6HPHI =»F13.7/1X»130 (1H-) 463 

5> '♦&'♦ 

2 F0RMAT(1X»*C0RRELATI0N(*»I3»*»*»I3»*)WAS OBTAINED USING A SUBSTAND '♦65 

2ARD CONVERGENCE VALUE OF*»F10,5) ^66 
RETURN 



467 



END '♦6^ 



SUBROUTINE APXTET3(PC»Pl'PJ»R0»EpS) ^69 

REAL M»K '♦70 

EPS=0. I'^l 
IF(PC.NE.PI)G0 TO 1 ' 

IF(PI.FQ.PJ)14»1 III 

TWOPI=&. 2831853072 ill 
M=PC-PI*PJ 
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IF(PI,GT.,5) PI=1.-PI 

IF(PJ.GT..5) PJ=1.-PJ '*77 

Ql=i.-PI '♦78 

QJ=1.-PJ '*79 

ZI=CDFNI(QI) '*80 

ZJ=CDFNI(QJ) '♦81 

H=EXP < . 5* ( -1 . * ( Z 1 **2 ) ) ) /SORT ( TWOpI ) '»82 

K=EXP ( . 5* ( -1 . * ( Z J**2 ) ) ) /SQRT ( TWOpI ) '♦83 

R=M/(H*K) '♦8*» 

R=R-(ZI*Zj)*(R**2/2.) '♦85 

IF(R.GT..80)U»2 '♦86 

2 lF(R.LT.-.80)3»5 '♦87 

3 R=-.80 $ GO TO 5 '♦88 
1* R=.dO '♦89 
"5 - RUT=ZI<^2-* ZJ«*2 - - - . - _s«90 

j=l f9l 

100 DO 11 1=1 » 50 '♦92 

IF(R.GE.1.)8»6 '♦g? 

6 IF(R.LE.-1.)7»9 '♦9^ 

7 R--.97 $ GO TO 9 '♦95 

8 R=.97 '♦96 

9 CONTINUE '♦97 
DIV=1.-R**2 '♦98 
FF'R I=EXP ( ( -1 . * ( RUT-2 . *Z I*Z J*R ) ) / ( 2 . *DI V ) ) / ( TWOPI *SQRT ( DI V ) ) U99 
F0FR=GAUSS6(R»ZI»ZJ) 500 
R1=R-(F0FR-M)/FPRI 501 
EPS=ABS(R1-R) 

lF(EPS.LT..0001)i5»lO 503 

10 R=R1 50*f 

11 CONTINUE ^05 
J=J+1 



506 



IF(J.6T.2)13»12 121 

12 R=.55 $ GO TO 100 5°® 

13 IF(EPS.LE.. 001)60 TO 15 =09 
R=99. $ GO TO 15 -^1° 

!*♦ Ill 

15 RO=R =12 

RETURN 513 

END l\'t 

FUNCTION 6AUSS8(R»ZI»ZJ) =lf 

DIMENSION U(a) »W(8) »U2(8) 516 
DATA (U(I)»I=1'8)/ .0198550718 , . lO 16667613 » .2372337950 » . U082826788 » 5l7 

+.5917173212'. 7627662050'. 8983332387*. 9801UU9282/ 518 
DATA(W(I).I=1»8)/.05061'*2681».1111905172».1568533229».1813'*18917» 5l9 

+ .1813'*18917'. 1568533229'. 1111905172*. 05061U2681/ 520 

DATA(U2(I) »I=1'8)/.00039U2239».010336130U». 0562798735* 521 

+.16669U7U58' .3501293883' .S81812283U' .8070026078' .96068U080'*/ 522 

SUM=0. if J 
TEMl=R/6. 2831853072 

TEM2=2.*ZI*ZJ*R 

R2=R**2 

RUT=ZI**2+ZJ**2 
DO 10 1=1 '8 

ADIV=1.-R2*U2(I) . c$n 

e=EXP ( ( -1 . * ( RUT-TtM2*U ( I ) ) ) / ( 2 . * ADI V ) ) /SQRT ( ADI V ) 530 



525 
526 
527 
528 
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10 SUM=SUM+6*wa) %li 
GAUSS8=SUM*TEM1 
RETURN 
END 
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